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Abstract 

Gauge-invariant non-local scalar and field strength operators have been argued 
to have significance, e.g., as a way to determine the behaviour of the screened 
static potential at large distances, as order parameters for confinement, as input 
parameters in models of confinement, and as gauge- invariant definitions of light 
constituent masses in bound state systems. We measure such "correlators" in the 3d 
pure SU(2) and SU(2)+Higgs models on the lattice. We extract the corresponding 
mass parameters and discuss their scaling and physical interpretation. We find that 
the finite part of the MS scheme mass measured from the field strength correlator 
is large, more than half the glueball mass. We also determine the non-perturbative 
contribution to the Debye mass in the 4d finite T SU(2) gauge theory with a 
method due to Arnold and Yaffe, finding (^m^i ~ 1.06(4)(7^T. 
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1 Introduction 



The strength and range of forces described by gauge theories are characterized by the po- 
tential of a static charge. Non-Abehan pure gauge theories in three and four dimensions 
are known to be confining in the sense that the potential in the fundamental representa- 
tion rises linearly with distance. If matter fields are added the static potential is modified: 
the linear rise persists only up to a critical distance at which enough energy is stored 
in the flux tube to pair create matter. The matter particles will then be separated and 
form bound states with the sources, so that the rise of the static potential saturates at 
the energy of the two bound state systems. If one considers static sources in the adjoint 
representation instead, the potential again rises linearly, but is supposed to saturate even 
in the pure gauge theory due to the creation of a pair of gluons. 

In this work, we investigate non-perturbatively operators related to the energy of 
the bound state systems created after the breaking of the flux tube in the pure SU(2) 
and SU(2)+Higgs models. The non-local gauge-invariant operators measured consist of 
sources of a given charge coupled via a Wilson line in the same representation: 



The fundamental Wilson line here is U^^^'^{x,y) = V exp{ig J^dxiA'^T"') and the adjoint 
Wilson line is t/„t^(x, y) = 2Tr T"^7^^"'i(a;, y)T\U^'''"^{x, y)]\ where = a''/2 and a" are 
the Pauli matrices. 

More specifically, the mass signal M extracted from these operators has a number of 
applications. First, related to the discussion above, M is naively expected to determine 
the asymptotic value to which the static potential saturates, ^(oo) = 2M (see below). 
The most prominent way of calculating the static potential is to extract it from the 
expectation value of large Wilson loops measured in lattice Monte Carlo simulations. 
In four dimensions (4d) this has been done for the fundamental representation ([0 and 
references therein) as well as for the adjoint representation (p[ and references therein). 
Analogous calculations in three dimensions (3d) may be found in [^0. However, no hard 
evidence for the saturation of the potential could be obtained in these calculations. One 
possible explanation is that the string in a lattice simulation does not break because of 
a potentially bad projection of the Wilson loop onto the hadronized final state. Further- 
more, Wilson loop calculations are very expensive in computer time. However, assuming 
the validity of the above picture of flux tube breaking, one may combine knowledge of 



G<f>{x,y) 
GF,ijki{x, y) 




(1) 



(2) 
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the static potential at small distances with its asymptotic value obtained from measuring 
M, to arrive at a rough estimate of the critical distance where the potential flattens off, 
i.e., the "screening length". 

Second, the asymptotic behaviour of G,j){x, y) constitutes an order parameter for the 
phase of the theory: in a Coulomb phase, 6*0(0, y) ~ exp(— m|?/|)/|?/|('^~^)/^, whereas in a 
confinement or Higgs phase, the behaviour should be purely exponential 

Third, the parameters related to the operator Gp^ijki{x,y) contain important non- 
perturbative information about the QCD ground state P] and may hence act as in- 
put parameters in models of confinement, such as the stochastic vacuum model H 
(for recent references see, e.g., [9-11]). This is of relevance also in 3d, since 3d theories 
can be regarded as laboratories for studying the qualitative features of confinement in 
QCD |l^, 0. The advantage of the 3d theories is that they are superrenormalizable and, 
consequently, exhibit very good scaling behaviour so that one can extrapolate results of 



lattice simulations to the continuum limit with an accuracy at the percent level [|14| , |I5 
One especially useful feature of the SU(2)-|-Higgs model is that the presence of the Higgs 
doublet allows a smooth interpolation between the non-perturbative confinement and the 
perturbative Higgs regimes |]16|, 0. The physical significance of the 3d theories stems 
from the fact that they constitute the high temperature effective theories of usual 4d 
theories in the framework of dimensional reduction [0, [1^ . 



Fourth, in analogy with the heavy - light quark system |]T9[, the mass parameters 
extracted from G^{x, y) and GF,ijki{x, y) can be viewed as determining the masses of the 
light dynamical constituents in meson-like bound states. This aspect can be of relevance, 
e.g., for the constituent models proposed to apply in the (symmetric) confinement phase 
of the 3d SU(2)+Higgs model p|, |2^. In this context, the Wilson line operators have 



also been conjectured ||2^ to explain the value of the propagator mass obtained from a 



simulation in a fixed Landau gauge [21 



Finally, related to the determination of the light constituent mass of a bound state, 
the operator Gp^ijuix, y) in the pure 3d SU(N) theory can be used to measure the lead- 
ing non-perturbative contribution to the finite temperature Debye mass in 4d SU(N) 
QCD [p2[. The connection to 4d physics is again via dimensional reduction. 



The purpose of this paper is to measure the operators in eqs. (|T]), (0) in 3d (in 4d, 
these operators have been studied on the lattice in [|^, ^). We improve upon the measure- 
ments of GF,ijki{x,y) performed for pure 3d SU(2) in [^, and extend the measurements 
to the SU(2)+Higgs case including now also the operator G^{x,y) (some properties of 
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G^{x,y) have been previously studied in [^). The discretization of the operators in 
eqs. (|l]), (0) is explained in Sec. 2. Some perturbative calculations are performed in 
Sec. 3. The technical details of the simulations are given in Sec. 4 and the numerical re- 
sults are presented in Sec. 5. In Sec. 6 we discuss the physical implications of our results 
for the various topics sketched above, and Sec. 7 comprises the conclusions. 



2 The lattice operators and the connection to the 
static potential 

The continuum Lagrangian of the 3d SU(2)+Higgs model is 



C = h^^F^ + {D,ip)^D,ip + mlip^ip + Uip^ipf. (3) 



On the lattice, one introduces = (v^ v'), where (p = i(72^* and a2 is the Pauli matrix. 
After a rescaling of 0, the lattice action is (Pjj is the plaquette) 

5[0,f/] = /3^^^[l-lTrP,,.]-/3^,^5:^Tr0t(a;)[/,(a;)0(x + z) 

X i<j ^ X i ^ 



In order to make contact with the calculations in [jT5| we fix the ratio of the continuum 
scalar and gauge couplings to be 



^ ^ ^ = ^ = 0.0239. (5) 



At this parameter value the Higgs model exhibits a strong first-order phase transition 
upon variation of jSn, or y = Tn\/ in continuum notation. We pick the same two 
points in parameter space as in namely y = 0.089 and y = —0.020, representing the 
confinement and Higgs phase, respectively. Once Pg is chosen, the appropriate values of 
Ph and (3r are determined by the "lines of constant physics" which govern the approach 
of the 3d theory to the continuum limit pS. Results for the pure gauge theory may be 



obtained by considering only gluonic operators and simulating at jSn = 0. 

Consider then the gauge- invariant two-point function in eq. ([l|). The path F between 
X and y could be for instance a rectangular one as shown in Fig. 1(a). However, it has 
been demonstrated by simulations in the 4d SU(2)-|-Higgs model that the coefficient 



M of the exponential decay of G^{x,y) is independent of R ||23[. Hence, we restrict 
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Figure 1: The operators measured on the lattice: (a) The correlator G tf,{x ^ y) with a funda- 
mental Wilson line (in practice we have R = 0, and the correlator is denoted by G^{T) ). 
(h) The correlator G p{T) with an adjoint Wilson line. 



attention to i? = in the following. Choosing T to be in the 3-direction the simplest 
choice representing eq. ([|) is then 

1 



-Tr 



<P\x)S{x,T)<P{y) 



with 



N-l 



S{x,T)= \{U^{x + n?,), y = x + N-?,, 



(6) 



(7) 



n=0 



and T = Na. In practice, we average this operator over the whole lattice in order to 
improve statistics. One can also choose fundamental charge operators other than the 0's 
at the ends of the Wilson line (see Sec. 4). The expectation (to be tested below) is that 
for large T, G^{T) should decay as exp (— MT) both in the Higgs and in the screened 
confinement phase [Q. 

The connection to the static potential and to screening is provided by employing the 
standard picture of string breaking through matter pair production. Assuming that the 
expectation value of the Wilson loop has a perimeter law behaviour for a large closed 
loop C, W{R,T) oc exp(-M|C|) = exp[-2M(i? + T)], the static potential, defined as 
V{R) = -limT^oo[lnl^(/2,T)]/T or 



V{R) 



lim In 

T-+00 



W{R,T) 
W{R,T-V 



obeys 



V{R 



oo 



2M. 



(9) 



Here M is assumed to be the same mass parameter as measured from eq. (^, when an 
optimal choice is made for the operator at the ends of the Wilson line. The physical 
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interpretation of eq. (^) is that at infinite separation of the static sources, the static 
potential consists just of the energy of the hadronized system with each static source 
binding a dynamical charge. 

Consider then the case of a static source in the adjoint representation. The Wilson 
line in the time direction now has to be taken in the adjoint representation, and the 
gluon field binding to it is described by a spatial plaquette. Hence we consider the 
gauge- invariant correlator shown in Fig. 1(b), 

Gf.MT) = (4Tr (P,,(x)T»)r^(r)Tr {PUy)T')) , (10) 

where 

T'^^T) = 2Tr (T"5(s, T)T^S^{x, T)) , (11) 
and 5* and y are as in eq. (0). This may be rewritten as 

GF,ijki{T) = (2Tr [P,,(a;)5P,Vi/)5t] -^TrP,,(x)TrP«(^/)) 

= (Tr [P,,{x)S {pUv) - Pu{y)) , (12) 

where the last equality holds for SU(2) only. In our simulation, we use Gi? in this last form, 
with the components ij = kl = 12; we omit the spatial indices of Gf{T) = Gf,i212{T) 
in the following. The plaquettes are replaced by the sum over all four spatial plaquettes 
of the same orientation sharing the end points of the Wilson lines, i.e. P{x) = Puix) + 
Puix — i) + Puix — 2) + Puix — i — 2). To improve on the statistics, we again average 
over the whole lattice. Adjoint charge operators other than Pij can be considered as well 
(see Sec. 4). 

The relation to the static potential may be taken over from the fundamental case. It 
is known that the static potential in this case also shows a linear rise due to flux tube 
formation between the static charges. As the flux tube now is in the adjoint representation 
it couples to the gauge fields and thus is expected to break even in the case of pure 
gauge theory. The corresponding final state consists of the static adjoint source binding 
a dynamical gluon, a system that has been termed "glue- lump" in the literature 0, |]. 



3 Perturbation theory 

Although the purpose of this paper is a non-perturbative measurement of the correlators 
in eqs. (|I]), (Q), let us in this section study these quantities perturbatively. The main 
motivation is to see how the mass parameters measured depend on the lattice spacing a. 
We will also make some other computations in the Higgs and confinement phases. 
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3.1 Scaling with the lattice spacing 



The mass parameters measured from eqs. (|T]), @ turn out to contain a divergent part 
oc g^ln{a). There is no counterterm in which to absorb this divergence so that, in fact, 
there is no meaningful continuum limit. Nevertheless, measurements with a finite lattice 
spacing can be useful for a number of applications, as we will see. Due to the fact that 
the gauge coupling is dimensionful in 3d, the divergent part can be determined with a 
1-loop computation. Moreover, it can be seen that the only such contribution comes from 
the Wilson line, depicted in Fig. ^ Following [^, let us thus consider this object. 



(a) (b) 

Figure 2: The Wilson line self-energy in the fundamental (a) and adjoint (b) represen- 
tation. 



Letting the Wilson line be an infinitely long straight line in the z-direction, expanding 
the path-ordered exponential to second order in the fields and taking the derivative with 
respect to z, one finds for the coefficient M in the exponential fall-off. 



M 



gl lim 



9l 



°o dz 

oo 



dy ' dx{A:ix,0)A':iy,0)) 



dx{Al{x, O)At(0,O)) (T^T" 



la/3 



d'^p± 



1 



(27r)2pi + m2,- 



(13) 



Here a = (3 are isospin indices not summed over, C2 = N for the adjoint representation, 
and C2 = {N^ — 1)/{2N) for the fundamental representation. The gauge field propagator 
was taken in a general gauge in the Higgs phase of the SU(2)+Higgs model with the 
mass my/- The result in eq. (0) applies also on the lattice when the integration range is 
restricted to (— vr/a, vr/a) and the momentum in the propagator is replaced with p]_ 
p\ = J2'^=ipf, where pi = (2/a) sin(api/2). The intergal in eq. (|T3|) can be computed and 
one gets 



n/a (pp 



-Va(27r)2p2 



32 



0{a) (lattice). 



(14) 
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/- — ,^ ^ — = — h In — - (MS in continuum), (15) 



(27r)2-2^ p-i + m 



which provides a relation between the two schemes. Eqs. (pTSD, (0) tell how the mass 
parameters measured on the lattice depend on the lattice spacing. 

In the following, we will want to consider quantities which do have a continuum 
limit. From eqs. (0), (|T^) , it can be seen that this is obtained by subtracting from M a 
divergent part: 



M'^M-gf^\J^\=M-gl'^\J'^-^\, (16) 



where Pg = 2A^/ {g^^^) for SU(N). Here the mass scale needed to define the logarithm in 
eq. ([T^) was chosen to be g1. The quantity M' is a continuum quantity in the sense that 
in the MS scheme with the scale Jt = g^, the exponential fall-off is determined by 

Mjj^ = M' + gl^-. (17) 
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6 6 6 6 6 i 

(a) (b) (c) 

Figure 3: The graphs contributing to the correlator G(fi{T) in the Higgs phase at 1-loop 
order. A dashed line is a Higgs field, a wiggly line a gauge field, a double line a ghost 
field and the cross denotes a counterterm. A filled circle is the shifted field (p, whereas 
an open circle is the quantum field if' . 



3.2 Higgs phase 

Let us then consider the correlator G^{T) in the Higgs phase of the SU(2)+Higgs model 
in some more detail. The Higgs field is thus written as <y9 = + <^ 7^ 0. At leading 
order, one gets just the constant 0^0. The graphs contributing at the next order are 
shown in Fig. 3. The other graphs are straightforward, so let us consider the Wilson line 
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contribution, (3.c). For arbitrary T, we get 

3 



<^0(^)l(3.c) 



fie^'mw, T) + fie^^mw, 0) + /(m^., T) - f{mw, 0) 

(18) 



2 



Here ^ is the gauge fixing parameter of an gauge and 

-,,j2kv^- '''' 
witfi tlie limiting values (S = 3.175911535625) 

V m 111/ p—i^AA 

fim,0) = ^-^ + O{am^), f{m,z)^'^^-'-—-. (20) 
Ana 47r 47t\z\ 

The gauge parameter dependent parts in eq. (0) cancel those form the other graphs, 
so that the resuh is exphcitly gauge independent. In the hmit T — > oo, eq. (pTsD is 
dominated by a hnear term, G^{T)\^^^^ Lp'^Lp^—MT), where M is from eqs. (0), (0) 
with C2 = 3/4. 

Including also the other graphs, the complete 1-loop answer for large T is 



T^ln 



IQir mwCL 



where mw,mH are the perturbative and Higgs masses, respectively. The 1/a-term 
on the 2nd row cancels the divergences from /(m//,0), f{myy,0), so that the sum is 
finite. At T = 0, on the other hand, the result is the first two rows plus a divergent term 
E/(27ra). 

In the Higgs phase for G^(T), one can thus see explicitly how a pure exponential 
behaviour arises, after the exponentiation into <^^<^exp(— MT) of the first row in eq. (pl|) 
(analytic arguments for the exponentiation can be obtained with the cumulant expansion, 
see, e.g., p6|). This is in contrast to the typical perturbative terms on the 3rd row 
containing pre-exponential corrections. It will be seen that the pure exponential is indeed 
what is observed on the lattice. 

For the vector correlator in the Higgs phase, the corresponding computation is much 
more tedious. At the continuum limit one may write 

P,,(x) ^ e'-'3^T-W (22) 
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M-<ij = ^ln^. (24) 



in eq. (|T0[) , and the leading term in the MS scheme is 

Gp(T) oc {F^MF^2{T)) ^ {N' - l)^e--^. (23) 

Taking into accout the next order, one gets a contribution exp(— M^'^-'T) from the Wilson 
line multiplying eq. (|23|), where, according to eqs. (p!3|), (p!^), 

27r amvi/ 

However, there are many other contributions as well (in pure 4d SU(N) in the continuum, 
the graphs have recently been computed in 0]). We will not evaluate these graphs 
explicitly. It will be seen that on the lattice the behaviour of Gf{T) is again favoured to 
be purely exponential rather than with a prefactor as in eq. (p3D, and that the exponent 
is ~ + M^'ij. 

3.3 Confinement phase 

Consider then the confinement phase of the SU(2)+Higgs model, or the pure SU(2) 
theory. From eqs. (|l5), (P^ it is immediately seen that a perturbative computation of 
the correlators is not possible: the tree-level mass parameter mw vanishes so that the 
IR-regime makes M divergent. One can only extract the coefficient of the logarithm 
which determines how the mass M scales when the lattice spacing a is varied. 

There is, however, the following meaningful computation one can make (These 
considerations are in complete analogy with those for the heavy - light quark system in 
4d, see |19|.) To be specific, consider the 3d SU(2)+adjoint Higgs theory. 



^ = i^l^l + ^(AAo)'^(AAo)'^ + ^{m'^yA^.A-, + ^A^(A«Ag)^ (25) 

where the mass parameter is large. Suppose one wants to compute the mass Mhl of 
the "heavy-light" bound state hi = eijkAQF^f^. A perturbative computation meets imme- 
diately with IR-problems |^2|, In addition, lattice simulations in the full theory are 



difficult since the requirement aMuL ^ 1 is very stringent, due to the heavy constituent 
mass m}^ (see, however, ||2^). On the other hand, perturbation theory does work for 



integrating out Aq, since there are no IR-problems due to the large mass. Therefore, 
one can integrate Aq out analytically, replacing {hi{0)hi(T)) by precisely the correlator 
Gp{T) in eq. (Q) in the pure 3d SU(2) theory, times an overall factor. The only subtlety 
is that the heavy mass Mh ~ appearing in the overall factor gets modified in the 
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integration procedure. As a result, the mass signal M measured from GpiT) represents 
the difference 

M = Mhl - Mh, (26) 
where Mh is determined by the computation in ||22| to be (for Aa ^ gi) 



^27) 

Then a lattice measurement of M from Gf{T) in the simpler pure SU(2) gauge theory 
allows a non-perturbative determination of Mhl, using eq. (|26|). These considerations 
are directly relevant for the numerical determination of the finite temperature Debye 
mass rriD in SU(N) QCD; we will return to this subject in Sec. 6.1. 

A computation similar to that for the adjoint Wilson line described above could be 
carried out for the fundamental Wilson line in the SU(2)+Higgs model, as well. This 
would allow a determination of the mass of the bound state if'^x in ^ 3d theory with two 
SU(2) Higgs doublets of which x is heavy. In the limit that (p and x interact only 
through gauge interactions, the computation is quite analogous to the one above, with 
the change (A^^ — 1)/(2A^) relevant for the fundamental representation. 

4 Simulations and analysis 

The algorithm used to perform the Monte Carlo simulation using the action in eq. (H) is 



the same as in |T^, [T^. The gauge variables are updated by a combination of heatbath 
and over-relaxation steps according to PD| , ^ , while the scalar degrees of freedom are 
updated combining heatbath and reflection steps as described in [^|. The ratio of the 



different updating steps is suitably tuned such as to minimize autocorrelations. In our 
simulations we typically gathered between 5000 and 10000 measurements taken after 
such combinations of updating sweeps. 

In order to investigate the scaling properties of the mass parameters and to extra- 
polate to the continuum we have simulated with Pg = 9, 12, 16 (/3g = 4:/gla). For each 
value of "we have checked that the extracted mass parameters are free from finite size 
effects, using lattices ranging from 14^ ■ 20 to 54^. 

4.1 Smearing and matrix correlators 

As mentioned in Sec. 2, the choice of operators at the ends of the Wilson line in eqs. (^), 
(p!oD is not unique. One could choose any other operators in the same representation as 
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well, and the task is to find those operators which give the smallest mass parameters 
(i.e., which couple to the "ground state"). A systematic way of doing this is to take a 
set of trial operators and to measure the whole correlation matrix, which can then be 
diagonalized. 

To form the basis of operators, one possibility is to take "smeared" or "blocked" fields. 
It has been demonstrated that using smeared variables instead of the original ones greatly 
improves the projection properties of gauge-invariant operators employed in calculations 



of the mass spectrum [|T5l. Similar findings have been reported from calculations of the 
adjoint Wilson line correlator in eq. ( |10|) P, ^. 



To be specific, we construct link variables of blocking level n according to ||32 



{x)ut'^\x + i) 



3 
+ 



±2 



E U^''-'''ix)Ur''ix + j)Ur''ix + l + j)Ur'^\x + 2i) ) ,(28) 



r(«-i), 



j{n-l) 



{«-l)t/ 



i=±i,jyi 



and composite scalar variables of blocking level n as in [|15 

2 



X] 



X 



(n-l)t 



{x - i)(t)^'''^\x - l) 



i=l 



(29) 

where i = 1,2, i.e., smearing is performed in the spatial plane. The Wilson lines con- 
necting the fields are in the time direction and remain unsmeared. Possible improvement 
techniques for the time-like links have been discussed, e.g., in 0, ^. 

The diagonalization of the correlation matrix is performed using a variational method 
(see, e.g., [l^, |l5l). For this purpose we iterate our spatial smearing procedure four 
times and measure the 5x5 correlation matrices 



GriT) = /^E^Tr [0Wt(a;)50(-)(y) 



(30) 

(31) 



where the plaquettes in P*^") have been constructed out of smeared links at blocking 
level n = 1, . . . , 5. For a given set of smeared scalar fields, we find the linear combination 
$1 = J2'^ik<P^''^ that maximizes G(j,{a), thus isolating the linear combination giving the 
lightest mass parameter. The first "excitation" may be found by repeating this step 
restricted to the subspace {0*^"^'} which is orthogonal to the ground state. Hence we can. 



11 



in principle, obtain five eigenstates of tlie matrix correlator given by 

<^i = j2a,k(p^'\ (32) 

k=l 

Our mass estimates are then obtained from the correlation functions (T) calculated 
in the diagonalized basis. Since the basis gets smaller for higher excitations, the reliability 
of the mass estimates rapidly deteriorates for higher states. Of course, one could improve 
on this by extending the basis of operators. However, even a small operator basis has been 
demonstrated to work quite well in the case of the lowest gauge-invariant eigenstates [1^ 



In the following we shall present results for the ground states and, where the reliability 
seems to be reasonable, the first excitation. 

Note that, in contrast to calculations of correlation functions of gauge- invariant oper- 
ators, the eigenvectors $4 themselves are not gauge-invariant and hence do not represent 
eigenvectors of the lattice Hamiltonian describing physical states. Nevertheless, the <l>i 
are gauge-covariant and the variational procedure will work to minimize the exponents 
of the fall-off. The same considerations hold for the field strength correlator (see also ) . 

4.2 Fitting functions and error analysis 

An important question concerns the fitting functions to be chosen in order to obtain 
mass estimates from the two-point functions. Unlike in the case of typical spectrum 
measurements where time slice correlators projecting on zero transverse momentum are 
considered, it is not a priori clear which asymptotic form the non-local operators take for 
large T. In particular, perturbation theory suggests that there are terms with exponential 



decay modified by power law corrections, see eqs. (|2ll ) , (|23|) . On the other hand, according 
to 0, one expects a pure exponential decay. In order to clarify this question we consider 
three fitting functions, exp(— MT)/T" with a as an open parameter, and the special cases 
exp(— MT)/T and exp(— MT). As we shall see, a pure exponential is the asymptotic form 
preferred by the diagonalized data (for the undiagonalized data the eventual conclusion 
is the same, but to achieve it requires a more tedious analysis). 

Our final mass estimates were therefore obtained by performing correlated fits of the 
form ~ exp(— MT) over some interval [Ti,T2] to the diagonalized correlators Gq,^^Fi{T). 
We have checked our results for stability under variations of the fitting interval, and also 
for compatibility with the results of uncorrelated fits. In those cases where different fitting 
procedures gave results that were not compatible within errors we quote the discrepancy 
systematic error. 
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exp(-MT)/T" 


exp(-MT)/r 


exp(-MT) 


Ti,T2 


aM 


a 


xVdof 


aM 


xVdof 


aM 


xVdof 


7, 12 
6, 12 
5, 12 


0.71(11) 

0.655(55) 

0.651(24) 


-0.30(84) 
0.17(35) 
0.16(14) 


0.63 
0.58 
0.63 


0.5366(96) 
0.5217(50) 
[huge error] 


1.15 
1.55 


0.6700(93) 
0.6758(46) 
0.6795(29) 


0.50 
0.51 
0.61 



Table 1: Comparison of fits to Gfi{T) in the confinement phase with different functions 
over various fitting ranges at (3g = 9, ■ Lt = 42'^. 



5 Numerical results 



In this section, we present our main numerical results. In Sec. p.l| we demonstrate the 
necessity of smearing and diagonalization in order to get a satisfactory signal with rea- 
sonable computer resources, and discuss the asymptotic form of the correlators. The 
mass estimates obtained for (3g = 9, 12, 16 and for various lattice volumes are contained 



in Sec. 5.2 and the extrapolation to the continumm limit is discussed in Sec. 5.3 



5.1 Asymptotic form of the correlators 

An example of the effects of the smearing and diagonalisation procedures is shown in 
Fig. ^ The field strength correlator is depicted in the confinement phase at /^c = 9 
on a 42^ lattice. This corresponds to the largest volume in physical units that we have 
considered. On the left the unblocked correlator is shown together with the once and 
twice blocked versions, all normalized to one at zero distance. It is immediately apparent 
that the unblocked correlator exhibits a faster decay at small distances than the blocked 
ones, and hence it is more difficult to extract the asymptotic mass value, although the 
data seem to be quite good. Furthermore, it turns out that all three correlators have 
a slight curvature. It requires a careful analysis of the stability of M with respect to 
different fitting ranges in order to decide that exp(— MT) gives a better fit than, for 
example, exp(— MT)/T. 

The situation can be significantly improved upon by performing the diagonalization 
outlined in the last section. The diagonalized correlation functions corresponding to 
the three lowest states are shown on the RHS in Fig. ^ In Table ^ we compare the 
results of fitting different functions to the diagonalized correlator with the lowest mass, 
Gp^iT). Acceptable fits using exp(— MT)/T may be obtained on the late timeslices. 
However, if one extends the fitting range to earlier time slices the quality of the fit 
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Figure 4: Le/i.' The correlator G f{T) for three different blocking levels in the confinement 
phase. Right: The diagonalized correlators GpiiT) for the three lowest states. Here Pg = 
9, volume = 42^ T = Na. 





exp(-MT)/T" 


exp(-MT) 




aM 


a 


xVdof 


aM 


xVdof 


Gpi, conf. 
G^i, conf. 
Gfi, Higgs 
Higgs 


0.655(55) 
0.3305(53) 
1.19(15) 
0.05424(8) 


0.17(35) 
-0.060(40) 
-0.4(5) 

0.0005(6) 


0.58 
0.53 
0.34 


0.6758(46) 
0.3228(13) 
1.09(2) 
0.05428(4) 


0.51 
0.72 
0.38 



Table 2: Comparison of fitting functions. In each case we consider the same fitting range 
for both functions. The data are for /3g = 9 from the largest lattices considered in each 
case. The last line corresponds to uncorrected fits and no x^/dof is quoted. 



rapidly deteriorates, giving large x^/dof and/or large errors. On the other hand, we 
obtain good fits with consistent mass values using a pure exponential exp(— MT) over 
a larger range of fitting intervals. Finally, three-parameter fits employing the functional 
form exp(— MT)/T" yield values for a that are consistent with zero as well as mass 
values consistent with those of pure exponential fits. The same behaviour is found for 
the scalar field correlator and in the Higgs phase as well. The fits for these cases are 
shown in Table ^. 

Thus, although there are a few cases where, e.g., a 1/T modification of the exponen- 
tial decay cannot be strictly ruled out based on a finite fitting interval and the value of 
X^/dof alone, c.f. Table |1|, we have strong numerical evidence that both of the correlators 
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G<s,(T) and Gf{T) decay with a pure exponential exp(— MT) for large T, in both phases. 
All the following mass estimates have been obtained by fitting a pure exponential to the 
diagonalized correlation functions G<i>^^F,(T)- 



5.2 Results for the masses 



The numerical results for the mass of the "ground state" in the Higgs phase are summa- 
rized in lattice units in the first block in Table ^. In it was found that the physical 
Higgs and W boson masses in the Higgs phase are free from finite volume effects already 
for rather small lattices, and our choice of lattice sizes is motivated by these findings. 
As a safeguard, we have performed an explicit check for finite size effects at /J^ = 9, 
and Table § shows them to be absent. The spatial volumes at = 12, 16 are larger 
than merely the scaled up versions of the smallest lattice for = 9, so we are confident 
to have reached the infinite volume limit in those cases as well. Note the very weak 
exponential decay of the scalar correlator. 
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0.401(4) 
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12 
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■42 






0.517(7) 


0.744(6) 




9 




30^ 






0.685(3)(7) 


0.960(10) 



Table 3: The coefficients of the exponential decay in lattice units. The first error is 
statistical and the second, where included, is an estimate of systematic effects. In the 
other cases the systematic errors were estimated to be smaller than the statistical ones. 



We remark that in the Higgs phase it was not possible to obtain any reliable infor- 
mation on excited states. In the case of the scalar correlator, the mass of the first excited 
state appears to be about an order of magnitude larger than that of the ground state. 
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However, the operators of all blocking levels have more than 95% projection onto the 
ground state, and correspondingly almost no overlap with the first excitation which is 
hence rather unreliable. In the case of the field strength correlator the effective masses 
computed from GpiiT) do not show a plateau, and correspondingly no fits are possible. 

The results for the mass in the confinement phase are shown in the second block 
in Table ^ Here, we have performed an explicit check for finite size effects for every 
value of Pg that we simulated. As is apparent from the table, at least all ground state 
masses have reached their infinite volume limits. At this point in parameter space we 
were also able to extract some estimates for the first excitations in the scalar and gauge 
field channels, which are labelled by Mq>^ and Mf2, respectively. 

Finally, the results for the two lowest states extracted from the field strength cor- 
relator in the pure SU(2) theory are given in the third block in Table ^ Note that the 
numerical values are very close to those found in the confinement phase of the Higgs 



model. This is not unexpected due to the by now familiar fact [|T5], |T3[ that the dynamics 
of the gauge degrees of freedom in the confinement phase of the Higgs model is rather 
insensitive to the presence of the scalar fields. The ground state mass for the pure gauge 
theory has been previously calculated in ||^. In comparing with that work we note that 
our result is ~ 5% lower. We ascribe this difference to a better projection achieved by 
employing the variational technique. 

5.3 Continuum limit 

Our next task is to examine the scaling behaviour of our results. The data are rewritten 
in continuum units in Table ^. The scaling of the data with jSc is shown in Figs. |^-|^. The 
observed slight increase of the mass values with decreasing lattice spacing is attributed 
to the logarithmic divergence discussed in Sec. [O]. Unfortunately, the divergence is so 
weak that in many cases it is not possible to clearly observe it numerically, based on our 
statistics and /^c-values. In any case, the logarithmic divergence has to be subtracted 
according to eq. ([16|) in order to obtain a finite continuum limit. The resulting mass 
values have been denoted with primes. For the primed quantities we observe rather good 
scaling consistent with linear ~ 0(a) corrections familiar from calculations of the physical 
particle spectrum of the theory |15|. (Let us note that for some observables the 0{a) 
errors could be removed analytically |Q, but such a computation has not been made for 



the mass values obtained from the composite operators in eqs. ([l|), (@).) The outcome 
of a linear extrapolation of the mass parameters in to /5g = cxd is given in Table ^ 
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Table 4: Masses in continuum units for the largest volumes in Table ^ (see also Figs. ^ 
1^. The primed variables are after the subtraction in eq. ( \r^ ) and are supposed to have 
a continuum limit. The linear continuum extrapolation and the corresponding x^/ dof are 
shown in the two rightmost columns. 

and constitutes our final result. 

In the Higgs phase one can now compare these results with perturbation theory. The 
tree-level perturbative W mass at this point is mw = 1-44(73. According to eqs. ([I6|) , 
(pID, the 1-loop perturbative scalar mass parameter is thus 

= -9l^ In - -0.022^3^ (33) 

which agrees within ~ 35% with the full value in Table |. It should be noted that a large 
logarithmically divergent 1-loop part has been subtracted in M^, so that this agreement 
is in fact quite good. For the vector mass M'p the order of magnitude expectation from 
the perturbative considerations in Sec. was M'p ~ mw Indeed we observe agreement 
in order of magnitude, but quantitatively there is a discrepancy. 
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Figure 5: The scaling behaviours of Mq>^, M$2 in continuum units. Left: Higgs phase. 
Right: confinement phase. The linear continuum extrapolations for the primed quantities 
are also shown. The logarithmic divergence of M^^ is so weak that it is not clearly visible 
in the confinement phase figure for the Pa-values available. 
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Figure 6: The scaling behaviours of Mp^, Mp^ in continuum units. Left: Higgs phase. 
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Figure 7: The scaling behaviour of Mp^ in the pure SU(2) theory. The linear continuum 
extrapolation for M'p^ is also shown. 

6 Discussion 

Let us now discuss the physical significance our results may have for the various topics 
outlined in the introduction. We begin with the application which has immediate physical 
meaning. 



6.1 The Debye mass 

An important concept in the phenomenology of high temperature QCD is the static 
electric screening mass, or the Debye mass m£,. However, it is a non-perturbative quantity 



beyond leading order |27|: the expression can be written as 



rriD 



m 



LO 
D 



+ 



In 



LO 



9l 



+ CNgl + 0{g'T), 



(34) 



where m]^ = {N/3 + Nj /GY^'^gT and Nj is the number of flavours. The logarithmic part 



of the 0{g'^) correction can be extracted perturbatively ||2^, but Cat and the higher terms 
are non-perturbative. To allow for a lattice determination, a non-perturbative deflnition 
was formulated in p2[, employing the dimensionally reduced effective theory of eq. (]25|), 
and the further reduction into the pure 3d SU(N) theory, discussed in Sec. 3.3. The 
statement is that the Debye mass can be determined from the exponential fall-off of the 
operator odd in which gives the lowest mass value in the theory of eq. (^); or, in 
terms of the pure SU(2) theory, from the exponential fall-off of an adjoint Wilson line 
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with appropriately chosen adjoint charged operators at the ends such that M obtains its 
lowest value. The latter is precisely the measurement we have made. 



In [22 1 it was further proposed that one could measure the constant M from the 
perimeter law of large adjoint charge Wilson loops instead of a single Wilson line, em- 
ploying essentially eq. (|^). As we have discussed, this measurement is in practice much 
more difficult than a straight Wilson line measurement: so far, it has not been possible 
to see the saturation of the static potential on the lattice in the adjoint case. 



Identifying now m^) with Mm in Sec. 3.3, it follows from eqs. 

Ar2 



cn 



Mf, N 
' H 



9l 



In 



1 



9l 



N 

+ —(In 4 



r 



71), (PI) that 



(35) 



Using the value M'pJ gl = 1.03(4) from Table |, we thus obtain ca = 1.06(4). 

The constant cat has previously been determined directly from the effective theory 
in eq. ( P^D for = 2, 3 in |^H[, where the physical consequences of the relatively large 
non-perturbative correction were discussed, as well. In [^], the result for = 2 was 
found to be C2 = 1.58(20). This is consistent with our result within ~ (2 . . . 3)cr, but our 
errorbars are much smaller. 



6.2 Constituent masses 

In ||2^ it was suggested that the non-local operators in eqs. (|l]), (^ would give a gauge- 



invariant handle on the masses of the light constituents forming the bound states in 
the confinement phase. The constituent masses are relevant for the bound state model 
of Ref . , as well. In it was further argued that the constituent masses should be 
consistent with the masses determined from the exponential fall-off of propagators in a 



fixed Landau gauge [^. This would then also correspond to the mass value obtained 



from gap equations in 



Let us first recall that within a constituent model, one would write the physical bound 
state mass Mhl as Mhl = Mh+Ml+Eb, where Mh and Ml are the masses of the static 
source and the dynamical charge, respectively, and Eb stands for the binding energy. 
Then in accordance with eq. (|26|), the mass parameter extracted from our correlators 
may be written as 

M = Mhl - Mh = Ml + Eb- (36) 

This still leaves the definitions of Ml and Eb open. In particular, it is important to realize 
that there is a divergence in M and thus in Ml or Eb- The precise way to handle this 
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Figure 8: The functional dependence on y at Pg = ^- Left: the lowest physical O"*""*" state 
(the Higgs mass in the Higgs phase) and the mass M' extracted from G$j(T). Right: the 
lowest physical 1 state (the W-boson in the Higgs phase) and the mass M" extracted 
from GpiiT). The central values of M" have been normalized by requiring M" = Mi — 
aty = -0.02. 

divergence cannot be decided without a further specification of the constituent modeL 

In it was suggested that a convenient way of circumventing the problem of the 
divergence and thus of defining a finite mass value M" (which should approximate Ml 
in some constituent models), would be to write M = M" + 5M with a constant 5M, and 
to require that in the Higgs phase, M" corresponds to the physical Higgs or W mass. 
This of course is only possible if the so defined M" has the same parametric dependence 
on, e.g., the continuum parameter y as the physical masses in the Higgs phase. If such a 
matching works, then M" will differ by a finite constant from the finite mass M' defined 
in eq. (|l^. 

The y-dependence of the mass parameters extracted from the field correlators is 
shown in Fig. ^ for /5g = 9. Clearly, the matching procedure does not work for the scalar 
correlator since the Higgs phase result does not behave as the Higgs mass with respect 
to y. The reason is that the signal in the Higgs phase is dominated by the cut-off effect 
found in perturbation theory, eq. (|21|). For the vector correlator, on the other hand, the 
parametric dependence of M" in the Higgs phase is found to be close to that of the W 
mass obtained from the 1 operator V^{x) ~ + ifj, and such a 

matching can be applied. Setting then 



M" ^M'j^^+Agl, 



(37) 



21 



A is determined from 



M"{y = -0.020) = mw{y = -0.020) = 1.91^^ (38) 

to be y4 = — 0.31(7|. The y-dependence of the M" thus obtained is shown in Fig.^ (for 
/5g = 9), together with that of Mi — . We observe that across the phase transition, M" 
jumps to a smaller value M"{y = 0.089) ~ 0-73g^ which then remains fairly constant, in 
contrast to Mi — . This functional behaviour is similar to that of the W propagator mass 



in Landau gauge as proposed in |20 



The numerical value of M" in the confinement phase may now be compared with 
those of other mass parameters extracted from lattice simulations. The lightest 0++ 
glueball mass is Mq = 1.60(4) (73 ||15|]; the fall-off of the W-propagator in Landau gauge in 
the confinement phase has been determined to be = 0.35{l)gl ||21]]. We thus observe 
that the numerical value of M" is about half the glueball mass and twice the propagator 
mass in Landau gauge. Thus, the conjecture in |Q according to which M" ~ 
appears not to be satisfied quantitatively. 

6.3 Screening 

The bound state system of a static source and a dynamical particle in the fundamental 
or adjoint representation is supposed to correspond to the hadronized final state of the 
respective static potential at infinite distance, see eq. (P). This permits, in principle, a 
rough estimate of the screening length Rg where the static potential flattens off, given the 
behaviour of the potential at small distances. The problem is that for the adjoint case, 
the potential has also been measured up to distances where the screening should show 
up if this picture is correct 0, but no sign of the screening has been seen. Nevertheless, 
it might be that the determination of the potential at large distances is less reliable than 
at small distances [Q. 

Assuming the potential indeed to be determined reliably only at small distances, one 
can obtain the following estimates from 

V{Rs) ^ 2M. (39) 

Since M is scheme dependent, one has to use consistently data from a fixed Pg- For the 
adjoint case in pure SU(2), our masses are ~ 5% smaller than measured in [^], but for an 
order of magnitude estimate this does not matter, and one gets ~ O-'^oi from Fig. 8 



in Q. For the fundamental case, we may use the fits made in ||, ^ to the Pg = 12 data 
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of i: 

V{R) =C- ^glKoimR) + aR, (40) 

where, corresponding to our confinement pliase point, the extrapolated fit parameters 
have been estimated as C = 0.310(71, m = l.05g^ and cr = 0.133(7| [^. Using the Pg = 12 
mass = O.TAOg^ from Table H, eq. (p9D then gives ~ 0.1(7| for the fundamental 
case. Given the uncertainties, these numbers are not to be trusted quantitatively, but 
it is nevertheless interesting to observe that both the adjoint and the fundamental case 
give a screening length of the same order of magnitude, and that the corresponding mass 
scale R~^ is quite small compared with all the correlator masses. 



7 Conclusions 

In this paper we have reported on the measurements of the expectation values of the 
operators in eqs. (|lD, (0) for Wilson lines with various lengths \x — y\, on 3d lattices 
with various lattice sizes and lattice spacings. This allows the extraction of mass signals 
related to the exponential fall-off of the expectation values, and a study of their scal- 
ing behaviour. Subtracting a perturbatively computable logarithmic divergence, we have 
been able to extrapolate the non-perturbative constant parts to the continuum limit. 
Apart from the mass signal in the exponential fall-off, we have also measured the asymp- 
totic functional forms of the correlation functions, verifying that they behave according 
to what is expected for confining theories. 

The mass thus measured from eq. (0) for the pure 3d SU(2) theory allows a deter- 
mination of the leading non-perturbative contribution to the finite temperature Debye 



mass niD with a method due to Arnold and Yaffe |Q. We obtain a much more precise 
estimate for rriD than has been achieved before: m£, = m]^ + g'^/{27^)ln{m]^ /g^) + 
1.06{4:)g'^ + 0{g^T) . The present measurement was made for SU(2) QCD. However, now 
that our study has proven the practical feasibility of this method, the extension to the 
realistic SU(3) case is straightforward. The physical significance of the relatively large 
non-perturbative correction has been discussed in 

On the more phenomenological side, the mass signals measured from eqs. (|l|), (0) 
are relevant for the composite models proposed to apply in the (symmetric) confinement 
phase of the SU(2)+Higgs model. This system has been studied as a laboratory for 
understanding confinement. The main observation based on our data is that, with a 
definite subtraction procedure, the lattice spacing independent part of the field strength 
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correlator mass is large: it is consistent with about half the glueball mass and twice the 
mass of the W propagator in Landau gauge. 

Finally, we have addressed the question of screening of the static potential. Assuming 
that the static potential measurements up to date are only reliable at small distances, one 
obtains both in the pure SU(2) theory and in the SU(2)+Higgs model the rough estimate 
1/Rs ~ 0.1 . . .0.2(^3 for the screening length of adjoint and fundamental representation 
charges. It is interesting that the mass scale here is much smaller than the other mass 
scales in the system. 
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